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Abstract 



The ARGO-YBJ (Astrophysical Radiation Ground-based Observatory at YangBa- 
Jing) experiment is designed for very high energy 7-astronomy and cosmic ray 
researches. Due to the full coverage of a large area (5600m^) with resistive plate 
chambers at a very high altitude (4300 m a.s.l.), the ARGO-YBJ detector is used 
to search for transient phenomena, such as Gamma-ray bursts (GRBs). Because 
the ARGO-YBJ detector has a large field of view (~2 sr) and is operated with a 
high duty cycle (>90%), it is well suited for GRB surveying and can be operated 
in searches for high energy GRBs following alarms set by satellite-borne observa- 
tions at lower energies. In this paper, the sensitivity of the ARGO-YBJ detector for 
GRB detection is estimated. Upper limits to fiuence with 99% confidence level for 
26 GRBs inside the field of view from June 2006 to January 2009 are set in the two 
energy ranges 10-100 GeV and 10 GeV-1 TeV. 

Key words: Gamma ray bursts. Gamma-ray, Extensive air showers, Cosmic rays 
PACS: 98.70.Rz, 95.85.Pw, 96.50.sd, 98.70.Sa 



1 Introduction 



Gamma-ray bursts (GRBs) are very strong gamma-ray photon emissions from 
cosmic unpredictable locations in a duration from milliseconds to tens of min- 
utes. They are the most energetic form of energy released from a single object 
in such a short time. The total amount of light emitted in a GRB is usually 
a factor of hundreds brighter than a typical supernova. Using thousands of 
GRBs detected by sateUite-based detectors, they have been thoroughly inves- 
tigated in the keV— MeV energy range. They are isotropically distributed in 
the sky with a non-thermal origin. According to the time duration, GRBs are 
usually classified into long (>2 s) and short (<2 s) bursts. Since the first de- 
tection by the BeppoSAX sateUite for GRB970228 [1] , afterglows are observed 
after GRBs are discovered, and this enables the multi-wavelength investiga- 
tion of GRBs from the optical band to X-rays. Redshift measurements show 
that GRBs occur at cosmological distances (the average redshift of GRBs ob- 
served by the Swift satellite is z = 2.3 [2]). Some short bursts come from inside 
old galaxies with little star formation, suggesting that they may be originated 
from mergers of binary neutron stars or black hole- neutron star systems [3]. 
Some long bursts are associated with supernovas and confirmed to be related 
to deaths of massive stars when central cores collapse to black holes [4] . 



* Corresponding author. Tel: -h86 10 88236106; Fax: +86 10 88233086 
Email address: chensz@ihep.ac.cn (S.Z. Chen ). 
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High energy (HE) gamma-ray emissions from GRBs is also observed in sev- 
eral satellite-born experiments. Energetic Gamma- Ray Experiment Telescope 
(EGRET) detected several GRBs with photon energies ranging from 100 MeV 
to 18 GeV [5]. Both prompt and delayed emissions were detected and no high 
energy cutoff was found in the spectra. Most importantly, a distinct HE spec- 
tral component was evidently detected in GRB941017 [6]. Recently, Fermi 
Large Area Telescope (LAT) also detected GeV emissions from GRB080916C 
[7] and 081024B [8]. These observations at high energies can place important 
constraints on models of emission processes and on parameters of the envi- 
ronment surrounding the sources of bursts. Very high energy (VHE) emission 
up to ~TeV is predicted by several models in both prompt and afterglow 
phases [9]. Emission at such high energies could result from electron Self- 
Synchrotron Compton (SSC) scattering in either internal or external, forward 
or reverse shocks. In such cases, a spectrum with double-peak shape extending 
into the VHE band is expected. Some models [10] also predict VHE emission 
due to decays of secondary 7r° mesons in neutron-rich outflows. Observations 
of VHE emission could play a role in discriminating between these models. 
The difficulty is that the absorption of the VHE photons by the Extragalactic 
Background Light (EBL), due to the pair production 7 + 7^51, —>■ e'^ + e~ , 
causes a substantial reduction of the VHE photon flux. This sets a high upper 
limit on the sensitivity of a detector used for GRB search in this energy range. 
The gamma-ray fluxes from these GRBs become too small to be detected from 
current satellite-based experiments due to their small sensitive areas, so only 
ground-based experiments have areas large enough for the detection. 

Search for VHE emission from GRBs has been done by many ground-based 
experiments including extensive air shower arrays and Gherenkov telescopes. 
No conclusive detection has been made up to now, while some positive indi- 
cations were reported. The Tibet ^^7 experiment found an indication of 10 
TeV emission in a stacked analysis of 57 bursts [11]. The Milagrito experi- 
ment reported evidence of emission above 650 GeV from GRB970417A with 
a chance probability of 1.5 x 10^^ [12]. Evidence of emission above 20 TeV 
from GRB920915G was reported about 1 min earlier than the GRB trigger 
time at 2.7a level by the HEGRA AIROBIGG array, but the position devi- 
ated of about 9° [13]. The muon detector GRAND found an excess during 
GRB971110 with a chance probability of 3 x 10"'' [14]. Due to hmited field 
of view (FOV), Gherenkov telescopes, like MAGIG and HESS, can only be 
operated in follow-up mode and at least 40 s (usually minutes) are needed to 
sway the telescopes to point to the GRB. However, this sets very low upper 
limits to the photon fluences at energies around hundreds of GeV during the 
afterglow phase [15,16]. 

With a large FOV (~2 sr) and high duty cycle (>90%), the ARGO-YBJ ex- 
periment, using a full coverage detector of Resistive Plate Ghambers (RPGs) 
with area 5600m^, is well suited for GRB surveying. Following alarms by 
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satellite-borne observations of GRBs, the ARGO-YBJ detector is used to look 
for emission from them with a threshold of a few hundreds of GeV. No sig- 
nificant excess has been observed yet. In this paper, we place upper limits on 
the VHE emission fluences for the GRBs inside the FOV of the ARGO-YBJ 
detector. Two models with very different high energy cutoff are used. The 
detector performance is investigated using MC simulation. Based on this, the 
sensitivity of the ARGO-YBJ detector for GRB search at different time du- 
rations and incident zenith angles is presented in section 3. The processes of 
the data analysis and the method to search for VHE emission are described 
in section 4. The results are reported and discussed in section 5. 



2 The ARGO-YBJ experiment 



The ARGO-YBJ experiment, a collaboration among Ghinese and Italian in- 
stitutions, is designed for VHE 7- astronomy and cosmic ray observations and 
located in Tibet, China at an altitude of 4300m a.s.l. The detector con- 
sists of a single layer of RPCs, operated in streamer mode, with a mod- 
ular configuration. The basic module is a cluster (5.7 x 7.6m^), composed 
of 12 RPGs (2.850 x 1.225m'^ each). The RPGs are equipped with pick-up 
strips (6.75 x 61.80cm^) and the fast-OR signal of 8 strips constitutes the 
logical pixel (named pad) for triggering and timing purposes. Hundred and 
thirty clusters are installed to form a carpet of 5600m^ with an active area of 
~ 93%. This central carpet is surrounded by 23 additional clusters ("guard 
ring" ) to improve the core location reconstruction. The total area of the array 
is 110 X lOOm^. More details about the detector and the RPG performance 
can be found elsewhere [17]. 

The RPG carpet is connected to two independent data acquisition systems, 
corresponding to the shower and scaler operation modes. With the scaler 
mode, each cluster of the ARGO-YBJ detector counts the rate of events that 
have total number of hits ^ 1, ^ 2, ^ 3 and ^ 4 every 0.5 s. Using the count 
rates, GRBs in the ARGO-YBJ detector FOV (zenith angle smaller than 45°) 
are investigated without direction information [18]. 

In shower mode, the ARGO-YBJ detector is operated by requiring at least 20 
particles within 420 ns on the entire carpet detector. The high granularity of 
the apparatus permits a detailed spatial-temporal reconstruction of the shower 
front. The temporal information is the arrival time of particles measured by 
Time to Digital Converters (TDGs) with a resolution of approximately 1.8 ns. 
This results in an angular resolution of 0.2° for showers with energy above 10 
TeV and 2.5° at approximately 100 GeV [19]. In order to calibrate the 18,360 
TDG channels, an off-line method [20] is developed using cosmic ray showers. 
The cahbration precision is 0.4 ns by using 24 h of data and the cahbration 
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result is updated every month [21]. 

The central 130 clusters began taking data in June 2006, and the "guard ring" 
was merged into the DAQ stream in November 2007. The trigger rate is ~3600 
Hz and the average duty cycle is higher than 90%. 



3 Simulation and sensitivity 

The effective area of the ARGO-YBJ experiment for detecting gamma-ray 

showers is estimated by using a full Monte Carlo simulation driven by COR- 
SIKA 6.502 [22] and GEANT3-bascd code ARGO-G [23]. Five zenith angles 
{6 = 0°, 10°, 20°, 30° and 40°) arc chosen in the simulation and the sampling 
area is 300 x 300m^ around the carpet center. The threshold multiplicity of hits 
(Nhit) which triggers ARGO-YBJ is 20. The effective areas for this threshold 
at the three zenith angles ^ = 0°, 20° and 40° are shown in Fig.l. 

10^ 
10' 

•s 1 

u 

£io'' 

10-= 
10 = 

10 10 
E,[GeV] 

Fig. 1. Effective areas of ARGO-YBJ for gamma rays with N^n > 20 as a function 
of the energy for the three zenith angles 9 = 0°, 20° and 40°. 

In order to avoid strong absorption of VHE photons by the EEL, two models 
of GRB emission with sharp cutoff of their spectra at 100 GeV and 1 TeV 
are investigated. Since the multiplicity of hits Nhu in an event is related to 
the shower primary energy, the optimization of the ARGO-YBJ sensitivity 
can be done by choosing a corresponding cutoff on N^it. As results of the 
optimization, ranges of Nhu are found corresponding to the two Ecut niodels 
as reported in Table 1. Given the ranges of Nh^, the corresponding optimal 
opening angle radii (pjQ, inside which 70% of the signal events is included, 
are reported in Table 1. This result is almost independent of zenith angle. 
After using the event selection listed in Table 1, the average effective areas 
(A) over the energy ranges from 10 GeV to 100 GeV and from 10 GeV to 1 
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TeV with a differential spectral index —2.0 (this assumption is used in all the 
following analyses) are calculated as functions of the zenith angle 9. They fit a 
functional form {A) = ylocos"6', where Aq = 4.36±0.04m2 and 119.3 iO.Sm^ 
are the average effective areas at ^ = 0°, with the parameter n — 14.26 ± 0.11 
and 10.56 ± 0.13 for the two energy ranges, respectively (see Fig. 2); with 
these functions the average effective areas at any zenith may be obtained. At 

9 = 20°, the effective area of the ARGO-YBJ detector is about 1.8m^ above 

10 GeV if the cutoff energy Ecut of the GRB spectrum is chosen as 100 GeV. 
On the other hand, it is 64:.3m^ if E^^t = 1000 GeV. 




5 10 15 20 25 30 35 
Zenith angle [degree] 



Fig. 2. The average effective area in two energy ranges after using the event cut 
listed in Table 1 as a function of zenith angles. The lines are the fitting result using 
the function cos" 9. 
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Fig. 3. The background rate detected in two different data selection conditions listed 
in Table 1 as a function of the incident zenith angle. This figure was made with data 
taken from one hour before to one hour after the GRB080903 trigger time, and every 
point stands for the average value from different azimuth angles. 



Using measured cosmic ray data, the background event rate (see Fig. 3) is es- 
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Table 1 

Nfiif ranges and corresponding angular window sizes for gamma rays. 



Ecut (GeV) 


Nhit 


070 (°) 


100 


20-60 


3.8 


1000 


20-500 


2.6 



timated according to the selected Nhu ranges reported in Table 1. Combining 
this with the effective areas shown in Fig.2, the minimum detectable fluences 
for GRBs requiring a 5a excess are estimated and shown in Fig.4. The sensi- 
tivity worsens with the increase of the zenith angle. The GRB time duration 
in this figure is fixed at T 12 s and the fluence scales with T^/^. We find 
that the ARGO-YBJ detector has a sensitivity of 10~^ erg/cm^ in the energy 
range [10,1000] GeV for GRBs with a duration of 12 s at ^ = 20°. 



10-^ 



10 = 



eio-* 



10-= 



lo- 



ll 0,1 00] GeV 
[10,1000] GeV 
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15 20 25 30 
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Fig. 4. The 5a minimum detectable fluences for GRBs with a duration 12 s in 
the two energy ranges [10,100] GeV and [10,1000] GeV for different incident zenith 
angles. 



4 Data analysis 



In an angular window surrounding the GRB with the aperture given in Table 

1, the Non events that fall in ranges are taken as on-source. Among 
these, the contribution from the cosmic ray background, N^, is estimated by 
integrating all the events over 2 h around the GRB trigger time, referred as 
"direct integral method" in the literature [24]. The significance of any excess 
in the on-source events is calculated using Eq. (17) of Li and Ma [25]. 

As widely discussed, the acceleration mechanism for VHE emission could be 
different from that at low energies, so that it could be on a different time scale. 
Even if the duration of every burst detected by satellite is known, the duration 
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Table 2 

The absorption factor K(z) due to the EBL is K(z) 
J A{E) {dN/dE)dE/J A{E) {dN/dE)e-^^BL ^E. 



FlnerP'v ( frpV^ 




^ = 0° 


9 = 10° 


9 = 20° 


9 = 30° 


9 = 40° 


10-1 on 

±.\J J-\J\J 


^ 


± . wo 


1 OR 

± . WO 


1 OR 

± . WO 


1 ns 

± .uo 


1 DQ 




1 


-L •KJtJ 


1 sq 


J- •KjrJ 




1 40 

J- .'rtKJ 




2.3 


3.39 


3.39 


3.40 


3.40 


3.54 


10-1000 


0.1 


1.60 


1.61 


1.66 


1.72 


1.75 




0.5 


7.92 


8.02 


9.27 


11.13 


11.66 




1.0 


25.71 


26.05 


31.59 


39.19 


41.84 




2.3 


116.5 


118.7 


147.6 


181.3 


209.2 



of VHE emission is still unknown. It is believed that (a) the most probable 
emission time is still in the prompt phase; (b) the emission may be delayed 
even by hours, like in the case of the 18 GeV photon in GRB940217, detected 
1.5 h after the prompt emission [26]; (c) the high energy emission could be 
produced earlier according to some models [27]. The GRB counterpart is first 
searched in a window Tgo, which is defined as the time in which 90% of the 
GRB photons is released. If no signal is found, we continue the search from 1 
h before to 1 h after the GRB. Time intervals of 1 s, 6 s, 12 s, 24 s, 48 s and 96 
s are used in this search with steps 1 s, 2 s, 3 s, 6 s,12 s and 24 s, respectively. 
If no signal is found, we set an upper limit to the fluence in Tqq. 

Given the observed Non and the expected background A*";, in Tgg, the upper 
limits for the number of events, Njjl, with confidence level of 99% is calculated 
with Helene's method (Eq. (10) in [28]). Using the effective area and assuming 
a differential power law photon spectrum, the upper limit to the fluence of a 
GRB in the VHE range is obtained. Guided by the average spectrum of the 
four bright bursts observed by EGRET, where a power law index of —1.95 ± 
0.25 was found over the energy range from 30 MeV to 10 GeV [5], an energy 
spectrum dN/dE = CE'"^ is assumed. At first, the normalization constant 
C is calculated by solving the equation Nul ^ ! A{E){dN / dE)e~'"'BL^E. 
Then the total fluences can be obtained by integrating / E{dN/dE)dE from 
10 GeV to 100 GeV and to 1 TeV. The tebl is the optical depth due to the 
EBL absorption, and the optical depths predicted by A. Franceschini [29] are 
used in this work. The absorption factor, which is defined as the ratio K{z) = 
J A{E){dN/dE)dE/J A{E){dN/dE)e-^^''^dE, for GRBs at redshift 0.1, 
0.5, 1.0 and 2.3 are fisted in Table 2. According to these results, the EBL effect 
for gamma-rays below 100 GeV is negligible, however, it could be substantial 
for farther GRBs in the range [10,1000] GeV. With the absorption factor K[z) 
and the average effective area (A) shown in Fig.2, the parameter C can be 
directly calculated by solving the equation Nul/{K{z){A)) = J {dN/dE)dE. 
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5 Results and discussion 



The data used in this work were collected by the ARGO-YBJ experiment in 
the periods July 2006-July 2007 and November 2007- January 2009. Twenty- 
six GRBs detected by satelhtes were within the FOV of the detector array 
while it was on. Angular resolution, pointing accuracy and stability of the 
ARGO-YBJ detector were thoroughly tested by measuring the Moon shadow 
in cosmic rays over all observational periods [30] . The Crab Nebula and flares 
of Mrk421 in years 2006 and 2008 were successfully detected [31]. 

No significant excess was observed for any of the 26 GRBs, neither as prompt 
nor prior/delayed emission (see Fig. 5). Upper limits to fluences in the VHE 
range are hsted in Table 3. The upper limits for the 6 GRBs with redshift 
information have been corrected for EBL absorption, while the others can be 
easily corrected using the factors listed in Table 2 assuming z= 0.1, 0.5, 1.0 
and 2.3. The fluences in keV bands measured by the satellite experiments are 
also hsted in Table 3. 



10^ F" 




Significance Significance 

Fig. 5. Left: Distribution of the statistical significance of the 26 GRBs with data 
available during their prompt phase using both events selections of Table 1. Right: 
Distribution of the statistical significance derived from two hours of observations 
around the 26 GRBs analysed with different time durations. The solid lines are 
normal Gaussian functions for given comparison. 

It is known that GRB spectra can be well fitted with a Band function with 
a break at energy E^) mostly between 100 keV and 1 MeV, and the average 
index a below the break is ~— 1 and the average index fi above the break is 
~— 2.3 [32]. Most GRBs observed by Swift do not have such a clear spectral 
structure with a break since its effective energy range is often lower than 
the break energy. GRB060805B was detected by the Inter Planetary Network 
(IPN) from 30 keV to 10 MeV and the result of fitting the data with the Band 
function [a = —0.66, f3 = —2.52 and Eq = 240 keV) is shown in Fig.6 (solid 
line). If this GRB spectrum extends to TeV energies only following the Band 
function, it will be compatible with the upper limits obtained by ARGO-YBJ. 
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Table 3 

List of GRBs in the FOV {9 < 45°) of ARGO-YBJ and 99% C.L. fluence upper 

limits 

GRB Satellite Redshift T9o(s) 9{°) keV fluence 10-100 GeV 10-1000 GeV 
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Fig. 6. The spectrum of GRB060805B and the upper limits obtained by ARGO-YBJ. 
The sohd Hnc is the result of fitting the IPN data with the Band function and the 
dotted line is a simple extrapolation. The upper limits are located at the mean 
energies of gamma ray events. For details about model 1 and model 2 see text. 

However, the spectrum might not extend with a simple power law (/? = —2.52) 
above 100 MeV. It might turn to (3 = —1.95 like the four bright bursts detected 
by EGRET observations [5]. If EBL absorption is negligible, the limit found 
by this work implies that the transition energy should be above 620 MeV (see 
model 1 in Fig. 6). The other possibility is that if the transition energy is 100 
MeV for this burst, the source should be farther than z = 0.15 (see model 
2 in Fig. 6). In addition, if the SSC mechanism used by Finke et al. ([33]) to 
interpret the spectrum of GRB 940217 also works on GRB060805B, the limit 
found by this work will provide a strong constraint. 

In conclusion, we have investigated 26 gamma-ray bursts in the field of view 
of ARGO-YBJ in about two years in the GeV— TeV energy range searching 
for prompt, delayed or prior emission. Since no significant excess is found, we 
have set upper limits to fluences for most of those bursts. 
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